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Abstract: The report is a first attempt to devise a calculation method 

for representing the buckling behavior of cylindrical shells 
of variable curvature. The problem occurs, for instance, in 
dimensioning wing noses, the stability behavior of which is 
decisively influenced by the variability of curvature. The 
calculation is made possible by simplifying the stability 
equations (permissible for the shell of small curvature) and 
by assuming that the curvature l/R as a function of the arc 
length s can be represented by a very few Fourier terms. 

We evaluated the formulas for the special case of an ellipse- 
like half oval with an axis ratio l/3 - e - 1 under compres- 
sion in longitudinal direction, shear, and a combination of 
shear and compression. However, the results can also be 
applied approximately to an unsymmetrical oval- shell segment 
under compression, shear, and bending so that the numerical 
values contained in the diagrams 10 to 12 represent directly 
dimensioning data for the wing nose. 
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1. INTRODUCTION 


So far, only the shell stability theory for sphere and circular 
cylinder has been developed to include formulas for practical application. 
Recently, the designs in airplane construction drew attention also to 
other shell types; however, due to the great mathematical difficulties 
opposing a stability theory of complicated shell forms, it was attempted 
here to approximate actual shells by circular cylinders. Thus it is 
possible, for instance, to calculate (with good approximation) a monocoque 
fuselage stiffened by stringers as if it were joined together from a 
large number of circular -cylindrical strips (shell segments). For the 
calculation of wing skins, one may use an approximation of the buckling 
formulas for the short circular-cylinder shell; the respective theory may 
be regarded as completed by Kromm’ investigations. Cases exist, 
however, where the theory of the circular cylinder does not suffice. A 
wing nose, for instance, is usually developed as a sort of half oval 
which is supported by a strong spar on the open side (fig. 1). The sta- 
bility behavior of such a shell is decisively affected by the variability t 

of the curvature. It is not permissible to determine the critical com- 
pressive load as in the case of the circular-cylindrical shell from a 
"mean value" of the curvature, since the weakening effect on the shell 
of regions of small curvature very considerably exceeds the stiffening 
effect of regions of larger curvature. (For rough calculation, it is 
even advisable to regard, in case of wide shells under compression, the 
minimum curvature alone as decisive for the buckling limit.) 

The mathematical difficulties opposing the theoretical investigation 
of cylindrical shells of variable curvature exceed those arising in the 
case of the shell of constant curvature; however, on the other hand, they 
are not of fundamental character for the infinitely long shell, the 
curvature variation of which over the arc length s may be represented 
by a sine series. The present investigation gives the general theory of 
such shells and evaluates the calculation results in particular for the 


%romm, A . : Knickfestigkeit gekrummter Plattenstreifen unter 

Schub- und Druckkraften (Resistance to Buckling of a Curved Plate Strip 
Under Shear and Compression) . Jahrbuch 19^0 der deutschen Luftfahrt- 
forschung, p. 832. 
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wing nose in the form of a symmetrical half oval. The following loads 
are considered: compression in axial direction of the cylindrical shell, 

shear (torsion), and the combination of both. 


2. THE STABILITY EQUATIONS OF THE CYLINDRICAL SHELL. THEIR 
SOLUTION IN THE SPECIAL CASE OF THE COMPRESSED 
SHELL OF CONSTANT CURVATURE 


The shell segment we are going to investigate is an intermediate 
between the flat plate and the true shell; it owes part of its load 
capacity to the reinforcement by the longitudinal-edge stiffenings (like 
the flat plate), part, however, also to the stiffening produced by the 
curvature (the curvature prevents a bending without stretching of the 
median surface). Thus the stability equations of the cylinder theory 
are buckling differential equations of the shell segment; however, their 
solutions have to be adapted to the boundary conditions at the edges 
s = constant. 


The two stability equations of the thin cylindrical shell are^ 


A ^ | Wxx = 0 


> ( 2 . 1 ) 


Et c 


lg ^ . v g - * - f Ax = “™xx + 2 J »xy 


where 0 is the stress function from which one obtains three "membrane" 
stresses additionally originating in buckling, according to 


Nx 

t 


< 3 2 0 _ 


«x ^ y _ j %s j Ng , 

— “ “ rsa> ~0xs> “ r xx > 

t 


a .is the critical compressive stress (hence the minus sign in equa- 
tion (2.1)), t the critical shearing stress; the remaining designations 
may be found in figure 2 which represents a shell segment. The indices 


2 Compare for instance K. Marguerre: Zur Theorie der gekrummten 

Platte. . ., Jahrbuch 1939 der deutschen Luftfahrtforschung, p. I^13. 
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appearing with 0 and w indicate the respective derivatives with 
respect to x or s. The equations each differ from the respective 

disk or plate equation by the j^-term. By precisely this term they are 

coupled with one another, and this is the reason for the mathematical 
difficulty of the shell problem. The equations (2.1) result from the 
very complicated "exact" shell equations3 by radical simplification.^ 
These simplifications are physically justified under the presupposition 
that the shell either has only very slight curvature5 (R » U) or at 
least is, in buckling, subdivided into so many waves that the wave length 
in the circumferential direction is small compared to the radius of 
curvature. If this presupposition is not satisfied, the result is to a 
certain extent problematical j however, the numerous careful investigations 
concerning the buckling of a plate of constant curvature^ have shown that 
the errors arising from the simplification of the equations are not of 
a type to endanger the technical usefulness of the results. 

In this report we consider only shells that are so long that the 
effect of end constraint on the buckling load can be neglected. Then 
the stresses and displacements are purely periodical in x, and for 

T = 0, in particular, directly proportional to sin — , with the wave 

length Z, which is small compared to the shell length L, remaining 
open at first. If we put 

0 = 0 sin w = w sin 

the partial equations (2.1) are transformed into total equations (bars 
now signify derivatives with respect to s) 



^Flugge: Statik und Dynamik der Schalen (Statics and dynamics of 

shells), Berlin 193^, P- 191. 

V few notes on this can be found in K. Marguerre's "Der Einfluss 
der Lagerungsbedingungen ..." (The influence of the arrangement condi- 
tions . . .). Jahrbuch 19^0 der deutschen Luftfahrtforschung, p. 867 . 

^A. Kromm: Die Stabilitatsgrenze eines gekrummten Plattenstreifens, 
(The stability limit of a curved plate strip), Luftf . -Forschg. 1938, 

P- 517 ff. 

^We name in particular: 

Flugge (footnote 3 ) 

Timoshenko: Theory of Elastic Stability, Chapter IX 

Kromm (footnote 1, p. 2) 

Donnell: Thin Shell Theory. Proc. of the 5th Inter. Congr. of 

Appl. Mech., Cambridge 1938, p. 66. 
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Et 2 

12(1 - v^V 





1 -T 

--o0 = O 

R Z 2 


( 2 . 1 ’) 


The equations (2.1') are linear and homogeneous in the unknowns w 
and 0 and contain the variable quantity l/R; they can be converted 
into a system of an infinite number of equations for the coefficients of 
a Fourier expression if l/R (or R) can be represented by a Fourier 
polynomial, and if the support at the longitudinal edges s = 0, s = U 
is such that by a Fourier expression the boundary conditions can, term by 
term, be satisfied. 


We consider first the cylinder of constant curvature because in this 
simple special case the train of thought leading to the solution is out- 
lined more clearly than in the general case; simultaneously, we have the 
opportunity of introducing the abbreviations which are expedient for the 
later calculation in an easily surveyable form. 

If we put 


w 


0 



4 ns 

an sin n -g- 



4 ns 

sin n -g- 



( 2 . 2 ) 


two equation systems linear and homogeneous in and A n are formed 

from (2.1') which with the abbreviation 


U/l = p 


(2.3) 


may be written 


^(p2 + n 2 ) 2 A n + | p 2 a n = 0 

(TT^y |( 02 + n2 ) 2a n - 


(2.4) 


= 0 


For l/R = constant = l/R 0 , only one n occurs in any equation, so 
that the system is broken down into pairs of equations only; calculating 
A n from the first equation and introducing it into the second, we obtain 
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^(p 2 + n 2 ) 2 _ 12(l - v 2 )u 2 ^ + 12(l - v 2 ) 


* 2 Et 2 


H 0 2 t 2 (p 2 + n 2 ) 2 ) n 


= 0 


The vanishing of the factor of a n characterizes the critical state. 
With the abbreviations 


0* - E xt 2 t 2 
1 - V2 3 xj2 


o> = 


12(1 - v 2 ) u 14 - 

X - ' 


R 0 2 t 2 


(2.5) 


the equation for buckling (by which the eigen-value a is determined) 
reads 

(P 2 + n 2 ) 2 


p^ , a 

+ <0 To p\ - ^ — = 0 

(p 2 + n 2 J a 


With the abbreviations 


1 _ 1 (p 2 + n 2 ) 2 

k n n/Hd P 2 




P 2 + 2n 2 + 


n^ 

P 2 




E t 


( 2 . 6 ) 


( 2 . 5 ') 


it is written in the still more compressed form 


— + k n - X = ° 


( 2 . 6 ') 


In equation (2.6) the wave length is still open in the x-direction 
^that is, p = yj. It adjusts itself so that a (or X) assumes the 
minimum value. Thus, there is, in addition to equation (2.6) or (2.6') 


da 

the condition = 0 or 

dX _ dX 
dp 2 ' d(l/k n ) 


(2 - 7) 
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From equations (2.6) and (2.7) it follows that 


P = n a = <x*|n2 + g j for \Z~ui ~ kn^ 


(e g + jl 

3 2 


= \/Zd a = i o*\faT 


(2.8)7 


or respectively K n = 1 X = 2 


? for \/7 d kn2 


The simple expression (2.2) is usuable only if at the longitudinal 
edges s = 0, s = U the equations 


w = 0, w ss = 0, N s = tj^ = 0, Eu = /(*„ - V0 XX ) dx = 0 (2.9) 

are prescribed as boundary conditions, that is, if the edge supports 
are developed in such a manner that they offer a very large resistance 
to a displacement of the sheet in radial and axial direction, a very 
small resistance to a displacement in hoop direction and to a torsion. 
For the further calculation, we shall assume the boundary condi- 
tions (2.9) to be satisfied. For the circular cylinder, a modification 
of the boundary condition for the stresses and displacements in longi- 
tudinal and hoop direction (thus relative to 0) can be shown to have 
almost no effect on the magnitude of critical load.® The same applies, 
with high probability, to the cylinder of variable curvature . For if- 
the edge terms are arranged in the more strongly curved region, they 
are practically without any influence because in this region the shell 
buckles into many small waves so that it does not matter at what point 
exactly the node lines are enforced} and if the stiffenings (as in the 
case of the wing nose) support the region of small curvature, it does 
become inportant whether the shell is simply supported or clamped 
(w- conditions) , but it is quite unimportant what happens in the two 


Tif one sets in the second o-formula the notations abbreviated in 
equation (2.5)* one obtains with v = 0.3 the well-known cylinder 

formula for buckling ~ 0.6Et 


a = 


R 


8 , 


K. Marguerre; compare figure 6 of footnote k, page k. 
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other directions, because for a flat shell (similar to a plate) large 
radial displacements w produce only very small hoop and axial displace- 
ments v and u. Moreover, the more ±mportant9 of the two assumptions, 
N s = 0 (bending soft supports in hoop direction), lies on the safe side. 


3. THE BUCKLING DETERMINANT OF THE COMPRESSED 
SHELL OF VARIABLE CURVATURE 


If the curvature l/R can be represented by a Fourier polynomial of 
r* terms-*-® 


- = —(l + 2^ cos — + 2{h cos — + 
R Rq \ U * U 




-U 


2r*ns\ 
2? 2r* cos —T'l 


u 


2rns \ 
COS — J 


(3.1) 


(because J cos ds = 0, l/Ro is therein the arithmetic mean 

of the curvature), we may, in order to solve equation (2.1* ), again 
start from the Fourier expression (2.2) which term by term satisfies 
the boundary conditions (2.9)* We obtain 


^Il( p2 + n2 f A n sln n f + | p2 H a n sinn 7 = 0 


U 


u 


, Et — p~r + n ?)% sin n ^ - cP 2 ^ sin n ^ V ( 3 .2) 

12(1 - V 2 )u 2 ^“ v ' u u 


Bin n f = 0 


^Only v (not u) is directly coupled with w over the hoop 

. . &v , w 

strain e v = — + „• 

y 8s R 

"*"^The equations of buckling for a similar case: i = — sin — 

R Ro U 

can be found in the yearbook report of the author quoted before. 
(Compare footnote 4, p. ^.) 
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whereby l/R ha8 to be substituted according to equation (3«l). On the 
basis of the trigonometric identity 


„ 2rits . nits . (n + 2r)its , . (n - 2r)its 

2 cos — — sm -jj- = sin ^ — + sin — 


(3-3) 


each of the two equations may be written in the form of a simple sine 
series: 


§ t(’ e * ^ ■*» =? ♦ f 


n=l 

sin ~ + g r)ns +JZa n sin 
r=l \n=l n=l 


U FL 


\ . nits 

Z_ a n sin — + 

n=l u 


(n - 2r)its 
U 


= 0, 


r (-^ — f — — ot(p 2 + n2 ) 2 - sin - ~ p £ 
STiM* 2 12(1 - / a r 0 


51 «n 


n=l 


^ njts , 
sin + 

U 


r* / 00 

Zl^2r(Xl A n sln 

r=l \n=l 


(n + 2r)ns 
U 


+ H An sln 

n=l 


(n - 2r)its 
U 


= 0 (3. If) 


By renaming the indices in the two last sums each equation can be made 

to contain only sin ~jy~* It is true that then the first sum begins 
at 1 + 2r ' and the last at 1 - 2r which involves a certain inconven- 
ience if one wants to read off in formally surveyable form the equations 
for the coefficients a n , A^ which result from the requirement that 
the factors of sin ££2. individually have to be zero. The difficulty 

can be easily avoided; one has 


r* / oo , v oo , x \ 

\ — f> K . (n + 2r) its . \ J (n - 2r)jrs\ 

Z_ ^2r 2— a n sin 1 tj- + A_ a n sln u" 

r=l \n=l X=1 / 


r* 


= TL Wll V 


-2r san “tt - + ^ — a n+2r 

r=l \ n=l+2r " u n=l-2r 


r 


sin 


nits 

U 


-A 


•t 
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or, splitting up the second sum with consideration of 
-sin (n - 2r)^ = sin (2r - n)^-. 



One can include the fully written terms in the first sum (which then 
starts with 1 instead of 1 + 2r), if one stipulates that by a k _ 2r 

(with negative index) the coefficient _a 2r-k (with positive index) 

must be understood, and puts a Q = 0. If one makes such a stipulation, 

all sums begin with 1, and (3**0 then reads simply 



Since the left sides of equations (3*5) for all values of the variables 
are to disappear, the curved braces each must be zero; thus one has 
the two equation systems for determination of the coefficients a^, A R 
of the expression (2.2) which, because of 

A -m = -V a -m = "V ( A o = a o = 0 ) (3.6) 

are valid for all n 

n = 1,2, 3 • • *oo 
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The first system does not contain the eigen-value cr, moreover 
only one A n appears in each equation; if we solve for A n we obtain 




l Is (B 2 ! n 2 ) 2 ( an + 


+ a. 


n+2r 


(3.7) 


Because of (3*6), this equation is valid also for n ^ 0; one can there- 
fore (without restriction for n) immediately substitute in the second 
system, and thus obtain one infinite equation system for determination 
of the unknown a n which, with the abbreviations (2.5')> can be 
written in the form 


(Vis - X H 


+ 5Z^2r( a n-2r + a n+2r]j + ^2 k n-2^ a n-2 + 


ZZ^2r( a n-2-2r + a n-2+2r)^ + k n+2^ a n+2 + ZZ ^2r( a n+2-2r + 


r* 


+ ^hj^n-U^n-^ + ^r(^-k-2r + a n-k+2r^ + 


1 n+2+2r ] 


( pj( V 

a n+U + 5Z ^2r( a n+4-2r + a n+i++2r ) j 


^6 [••] + •• • - 0 

(3-8) 


The determinant A of the homogeneous system (3-8) is the buckling 
determinant; A = 0 is the qualifying equation for the eigen- 


value X = 


a* \fw> 


a . 
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Equation (3-8) is somewhat modified by writing out the sums over 
r and rearranging them 

+ ^[( k n + V2H-2 + ( + ^-2)^2] + 

C 2 2 [ k n-2V-4 + k ,A4 + ( k n-2 + k n +2 K] + ? 4 [K + k n-4) a n-4 + 
( k n + k n4)v] + E 2 ? l.[ k n-2( a n-6 + Va) + k n-4( S n-6 + a n-2 ) + 
k n + s( a n-2 + Vfi) + k n + 4(\ + 2 + Vfi)]* V^nVn-S + k n + 4 a „.8 + 
( k n-4 + \ + 4) a n] + ? 6 [k + \sK-6 + ( k n + N^Vg] + 

?2?6 k n- 2 ( a n-8 + a n+4) + k n-fi(V8 + Vu) + + V-fl) + 

k n+6 ( a n+4 + a n +8 )j + ?4^6 [V-4 ( a n-10 + a n+ 2 ) + k n-6( a n-10 * a n-2 ) * 
k n+4( a n-2 + a n+lo) + k n+6( a n+2 + a n+lo)] + j^n-6 a n-12 + 

k n+ 6 a n +12 + ( k n -6 + k n-t- 6 ) a nJ + • • • * 0 ( 3 . 9 ) 

The system (3-9) contains only even or odd indices - that is, the 
buckling forms symmetrical and antimetrical with respect to the center 
do not effect one another (of the two buckling forms, the one to which 
the smaller buckling load X pertains will occur) . 

If the series (3*1) breaks off, as we assumed, after r* terms, 
we are dealing with a finite-term system. If (3.1) contains, for 
instance, only the two first terms, a n _^, a n _ 2 , a n , a^g, a n+ 4 > 
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appear in every equation - the system has five terms. For ^ 0 
it becomes a 9-term system, for f 0 a 13-term system, etc. 

One determines X by first completely setting up (as was done 
here for r* = 3) the system for the given number r* of £ -terms, 
then arranging it according to the unknowns a n (for instance 

a n-2^2( k n + k n-2) + %( k n+4 + k n-6)J etc-) and calculatin S the 
coefficients of the a n which thus originate. For prescribed £ -values 
these coefficients are, furthermore, functions of k n , that is of cn, 

P, and (in the main diagonal) of X. The curvature value cd is to be 
regarded as prescribed, (3 is chosen, and X is then determined from 
the condition that the coefficient determinant must disappear. The 
suitable 3 would have to be determined - as in the special case 

l/R = const. - from the condition ^ = 0 which expresses that the 

wave length l = U/3 appears pertaining to the smallest buckling 
value X; however, due to the high degree of the X-equation, this 
condition can be satisfied only by plotting X(P) for a number of 
neighboring P -values and reading off an approximate X TTli r| . 

The method requires much detailed calculation - still, it is 
superior to other methods (for instance to the method of reducing the 
number of equations by a to-be-guessed relation among the coefficients 
in the expression (3*1) - which can be interpreted as a sort of Ritz 
method), because it can be highly systematized and, above all, because 
its accuracy may - starting from rough approximation values - be 
increased arbitrarily and at any time. 

Nevertheless, the calculation iB rather troublesome, due to the 
large number of parameter values co, (( for which it must be performed. 
In order to make a useful choice among the many possible parameter 
combinations we must therefore, before numerical evaluation, answer the 
question as to the appearance of the shell forms corresponding to a 
certain variation of curvature l/R(s). 

k. DETERMIMTION OF THE CYLINDRICAL FORM 
FROM THE VARIATION OF CURVATURE 


By the equation (3-1) 

1 1 L „ 2ns , n y 4ns y 6ns 

R = R^\ + 2t >2 cos — + 2C 4 C0S IT + 2>i 6 COS ~U 


• .) (4.1) 
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the curvature is given as a function of the arc length. Because of 

1 _ dffl 
R ds 

it follows therewith that for the angle cp = cp(s) 

P ds s „ U 2 its „ U . 4its „ U . 6ns , 

cp = / — = — + Co -4- Sin + Ck — 1 — Sin + Cg sin + • % 

J R R 0 JtR 0 U 2nR Q U ° 3*R 0 U 

( 4 . 2 ) 


From (4.2) one obtains by numerical integration because of the relation 


dy dx „ 

- 7 *- = sin cp, — = cos cp 

ds ds T 

the ordinate y and abscissa x of the profile curve: 


(4.3) 


y = 


sin cp ds, x 


cos cp ds 


(4.3') 


The ratio U/R q in equation (4.1) is at first still open. It is 
determined from an assumption regarding the opening angle cp Q = cp(U) - cp(0) 
of the shell; according to equation (4.2) it is 


<Po = 



We shall study here, above all, the half oval. (Compare fig. 1.) 
For it 


cp Q = it, that is, U = RqJT 
and equation (4.2) is transformed into 


(4.4) 


SIT 


2its 


4jts 


cp = — + Co sin ~~ + — Cli sin 
U d U 2 4 U 


+ — C/T sin 


6its 

u 


+ 


(4.4') 


In order to reduce calculation expenditure in application of the 
equation system (3-9)> it is desirable to manage with as few terms in 
the series (4.1) as possible. In figures 3a and 34 curvature and course 
of curves according to equation (4.1) from (4. 4') are plotted for the 
parameter values Cl^. = Cg = ••• =0 and C2 = -0.6 or -1_, respectively. 
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Figure b contains the curves required for conversion from the parameter 

b b b 

£ = -£p to the axis ratio e = — and to the quantity ? ' = — = — . 

a R o U 

It is seen that one can just reach the axis ratio — = i with the two- 

11 a 2 

term expression: for values |^ 2 |> — one voul(i have a change of sign 

in the curvature course, that is, one would obtain an oval buckled 
inward at the ends of the small axis - a form that does not occur in 
the applications. 


In order to obtain "reasonable" ovals with an axis ratio € < 1 / 2 , 
one must therefore start from a multiterm expression of the type (4.1) . 
Figures 5a and 5 h show (as the most obvious oval type) an ellipse of 


the axis ratio 



together with its curvature variation. 


The 


figure explains why the simple expression 


i = -Lfl - i cos 2*s\ 
R R 0 V U / 


(U.5) 


in the region of 



that is. 


s = 1 may no longer produce ellipse- 


like curves. Simultaneously, it indicates what type of curvature 
expression must be chosen: For the limiting curve (the oval with 

vanishing curvature at the end of the small axis) a higher power of 

JtS 

sin — comes into consideration for 


?6 = ?8 - 


. = 0 , thus, for instance. 


I ~ sin 1 * ££ 
R U 


that is 


For £3 = £ 10 


1 

R 



b 2 ns 1 
3 cos y + 3 


cos 



. . . = 0 


1 

R 




sin 


6 ££ 
U 


( b.6) 


1 

R 



15 

10 


cos 


2 ns 6 
U 10 


cos 


4ns 

U 


10 


cos 



(4.7) 


that is 
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Thus one will use for the pertaining oval -f ami lies for instance the 
expression 


thus 



s(l -33 COS - 0.33 COS 


^2 = -1-335, £4 = +0.335, 


or, respectively. 


1 

R 



5^1.5 cos — - 0.6 cos 


Uns 

U 


+ 0.1 cos 


6irs\ 

U / ' 


(4.8) 


( 4 - 9 ) 


thus 

5 2 - - 1 . 55 , t k = + 0 .65, 5 g = -0.15 


Figure 6 represents the oval 5 2 = _1 from figure 3h together with the 

oval determined by equation (4.8) with 5 = 0.8 and the ellipse of the 
same axis ratio. It can be seen that the oval corresponding to the 
multiterm expression is more ellipse-like, in particular, that it shows 

no position of zero curvature. In figure 7 > e = — , e' = — , 

a U 


and 


P = 


3/R(uZ2) 

r(o) 


are plotted for the expression (4.9) • 


For an ellipse 


p would be equal to € - one recognizes from the representation that 

p( 5> and e(£) between e = 1 and € = 0.35 lie close together so that 
the four-term expression (4-9) is sufficient to describe reasonable 


ovals down to an axis ratio 



1 

3 ’ 


In the expressions (4.8) and (4.9) the signs alternate. If one 
chooses the same expressions all with the same signs, it signifies that 

JT 

the positions cp = 0 and 9 = 2 exchange their roles: One obtains 

(with exchanged axes) the same oval; only now the two halves are 
connected to the flanges at the ends of the large axis, if we still 
visualize the flanges at cp = 0 and n. Since this case is of less 
practical interest, we shall not discuss it more closely - but it should 
be pointed out at least that it is in this simple manner related to the 
’’main case” equations (4.8) and (4.9), respectively. 
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For reasons of calculation expenditure we refrained in the present 
investigation from extending the expression for l/R stili further; 
going still further into an individual case does not present any funda- 
mental difficulties; however, the expenditure will hardly ever be worth 
while for the applications since the calculations performed with the 
expressions (4.5) to (4.9) led to the result that the axis ratio 6 
is a suitable parameter for characterization of the stability behavior 
of an oval, that is, that the form variation in detail is of minor effect 
on the magnitude of the buckling load. 


5. COMPUTATIONAL EVALUATION OF THE BUCKLING 
DETERMINANT IN THE CASE OF COMPRESSION 


Calculation of the eigen-values X from the condition that the 
determinant of the homogeneous equation system ( 3 - 9 1 ) must vanish is 
possible only "step by §tep:" Step by step in the sense that one starts 
with small t, and o> values and with the two-term expression (4.5) 
for the curvature variation. 


For = ... 0, each of the equations (3*9) contains 

only the two first lines, and the system reads, arranged according 
to the unknowns (with = “£)> 

for odd n 



kp(l + K ) 2 + £^3 - X 


a l 


£ (l + £)ki + k^ a^ + 


= 0 
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and for even n 

+ k 2 + ? 2 

-?[ k 2 + \] a 2 
?2 V 2 - ? [> 

^ ?k 6 a l* - 5 [kg 

The second system shows the usual initial irregularities in the first 
equation, the first system, those in the two first equations. For the 
first system these irregularities are such that the sign of £ 
considerably affects the behavior of the first terms. 

In contrast, the eigen-value resulting from the second system does 
not depend on the sign of £ (only a^, ag, a-^ change their signs); 

since the even terms in the expression (2.2) have nodes at s = 0 
and u/2 (cp = 0, jt/ 2), it must indeed be a matter of indifference at 

which of the axis ends the flanges are connected. 

The secondary terms are small if £ is small and - since they 
contain only the quantities ~ /nT - if /oT is small. For 

sufficiently small values of £ and /uT (for instance, 2^ = 0.5, 
y/co = 50), therefore, a portion of two equations from one of the two 
systems (5*1 or 5-2) suffices for the determination of the eigen-value. 

The quadratic equation for X is easily solvable, the minimum value 
as a function of P can be determined immediately. If one retains 
for instance £ and increases one can at first still manage with 

two or three equations, can determine p and X, and plot both 
quantities as a function of \px>. By extrapolation one obtains approxi- 
"mation values of P and X for the next /oo value; one improves these 
approximation values by trial substitutions in the determinant, and 
continues in this manner up to the upper limit for /n>, at about 
\Jla = 4000, which is of practical interest. 

A twofold difficulty opposes performance of this calculation: 

First one must decide on the order of the determinant portion to 

be selected which fixes with sufficient approximation the eigen-value |i 
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(which, actually should be determined from an infinite determinant Aoo) . 
Since the order of ^ increases with (the secondary terms become 

more and more important), the buckling determinant finally becomes of 
too high an order to be still solvable with respect to X. Neglecting 
the convergence problem (regarding the mechanically reasonable result 
as a sufficient confirmation of convergence), we can surmount both 
difficulties by the following method. Besides the one for 3, we 
insert an approximation value Xp - (estimated by extrapolation) for 

X also, and reduce the determinant to its main diagonal with the 
aid of Gauss' algorism - without first fixing m. If we plot the new 
diagonal values (the product of which represents the value of the 
determinant against the pertaining number of equations n, we 

obtain a point sequence we can connect by a curve. If we continue 
sufficiently far with m, this curve will either intersect the axis 
(generally not precisely at an integral n) or bend up. If the curve 
intersects the axis, the main diagonal terms were too unimp ortant, 
thus X = Xp too large. The calculation is repeated with a slightly 

smaller X = Xg) if the curve this time bends up, the number m was 

chosen correctly, and the desired X-value must be between Xp and X^ 

because the new diagonal terms would, with increasing n, have to 
decrease asymptotically toward zero (a» — > 0) for this X-value. 

Figure 8 shows a number of such curves for 2 (; = 0.5 and several sets 
of values for \/(E, 3, and X. From the two curves corresponding to 

\/~a> = 500, 3 = l6, one recognizes, for instance, that the root of A» 

must be between Xp = -l.l8 and X 2 = 1.16. The difference between 

the values Xp and X 2 is a measure of the accuracy of the calculation. 

The entire procedure must now be repeated for neighboring 3-values in 
order to make reading off the minimum from a X(3) curve possible. 

Since, however, X as an extreme value in the neighborhood of the 
correct value is affected very little by the choice of 3* two - at the 
most three - 3 proofs will generally be sufficient, provided the 

investigator has had some "experience" with this calculation method. 

One thus obtains a family of curves X ( \Zls) with £ as parameter. 
Since the curves run smoothly, one can conveniently interpolate on' the 
basis of four such curves (2£ = 0, 0.5> 0.8, 1.0) and with the aid of 
the curves of figure 4 convert to the directly prescribed quantities, 
the axes a, b of the oval, and to the wall thickness t contained 
in cr* and /oT. According to (2.5) 


a _ -v a* v/oT. 

— - A, — — ) 

E 4E 
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according to (4.4) one obtains for the half oval 




, / — 2\/ 3(1 - v 2 ) U J . 1 a X 1 Jtt t c -*\ 

V£D = — ^ > thus — = — — = X — r»; - — (5-3) 

ic t E V cd 2\/3{l - v 2 ) U 

tcId 

With the dimensionless e' = — represented in figure 4 one may write 
instead 


a 

E 


1 

1/3(1 - V 2 ) 



(5-4) 


and since 
(that is, 

the ratio 


X is itself a function of /cd (that 
U/b), everything can be converted to 

— (spar height/skin thickness), and 


of the oval. 


is, 

the 

the 


U/t) and £ 
two parameters, 

axis ratio — = 
a 


e 


In the region 200 ^ V<d <; 4000 which is interesting in practice 
(that is, for the wing nose), for ((-values >0, the two 1- values are 
not distinguished according to equations (5*1) and (5«2). This becomes 
physically understandable if one considers figure 9 where for a "mean" 
pair of parameters /cd, £ the buckling profile, as it results from 
equation (5*2), is plotted against the developed width. One can see 
that a large wave originates in the flat region while many tiny waves 
develop in the region of large curvature; and it is immediately plausible 
that it cannot affect the magnitude of buckling load whether these 
small waves happen to have a bulge or a node at the point cp = H. (For 

2 

values < 0, of course, the buckling values differ considerably; for, 
whether or not a node is enforced at the flat place cp = 0 is decisive 
for the buckling form. ) 


The calculation with the multiparameter curvature expressions (4.8) 
and (4.9) takes fundamentally exactly the same course, except for the 
fact that here the computation of the coefficients of the equation 
system (3-8) as functions of £ and of the quantities is much more 

troublesome, and that the complexity of the system is of disadvantage 
also for application of Gauss' algorism. With a four-term 
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expression (4.1) as a basis, the first coefficients of the system read 
(we denote them here simply by their indices put in parentheses) : 


(id - i ♦ *,(1 - y 2 ♦ * 3 (e 2 - e 4 ) 2 * * 5 (e 4 - e 6 ) 2 + y 6 ‘ 


(22) - 4 + M 1 ■ 4) 2 + ' ? s) 2 + k s^ 2 + k sV 


<33) = cj + ^(1 - y 2 * h(«2 - y 2 + k 5 ? 2 2 + n? 4 2 + ye 2 
(44) = ^ + \ + k 2 ( ? 2 - kf + V/ + k 8 ? l. 2 + k 10 ? 6 2 

(i3) - k i( ? 2 - M 1 - h) + k 3 ( ? 2- yt 1 - y + v 2 (v y + *Mh 5) 

(24) = k 2 (t, - y(i - e 4 ) + ^ - Eg) + - kgE^ 

(35) = ^(E, - E 4 )(E 4 - Eg) + k 3 E 2 (l - Eg) * *5E 2 + y 2 E 4 + *^6 
(15) = k 1 (E 4 - Eg)(i - E 2 ) - k 3 E 2 (E 2 - E 4 ) + k 5 ( ?4 - Eg) + t^Eg 


Fortunately, the calculation with the new coefficients need not be 
performed for all former \f(E and £ values j the (3 -values and 
the much less sensitive X-values for the ovals of equal axis ratio 

and equal ratio ~ are found to lie very close together. Thereby the 

calculation with the expressions (4.8) and (4.9) assumes the character 
of a check calculation for the region 1 ^ € ^l/2 and of an extra- 
polation calculation for the values of l/2 ^ € > 1 / 3 . 


That the buckling load depends only on 2b/t and b/a (that is, 
on the form "on the whole"), but not on the curvature variation in 
detail, can be physically explained by the effects of curvature varia- 
bility which oppose one another: The parts of stronger curvature support 

those of small curvature, and - as the calculation shows - it does not 
make much difference for the load capacity whether a small region with 


22 


NACA TM 1302 


strong curvature is present, which for the rest has the character of an 
almost rigid edge support, or whether, as for the oval (it-. 5 ) } a larger 
region of medium curvature expands which relieves the rest of some of 
its load. However, it is noteworthy that not the mean curvature (one 
parameter cn) is decisive for the critical compressive stress, but that 
two parameters, axis ratio and spar height, appear as the critical 
quantities . 

The total result is shown in figure 10 in logarithmic representation. 
Use of the curve table requires a brief experimental method since for 
reasons of clarity in the extensive U/t range the expression propor- 
tional to the compressive force 


£ JL 

E 2b 


(5-6) 


could not be represented, but only the quantity proportional to the 

stress Thus one must, for prescribed compression force, first 

estimate the quantities a or t and may then determine on the basis 
of the curve in what direction the estimated value must be changed. 


6. THE CLOSED CYLINDER UNDER PURE SHEAR 


The shear problem differs considerably from the compression problem 
in mathematical respect. The main equations (2.1) no longer contain in 
all terms an even number of derivatives with respect to each of the two 

TTX 

variables. Consequently, sin — can no longer be factored so that 

the equations (2.1'), which served as starting point for the compression 
calculation, lose their validity; one must refer back directly to the 
partial differential equations 


A A0 " 1 w xx = 0 


Et 2 

1 - V2 


A Aw + 




XX 


2tw. 


xs 


= 0 


( 6 . 1 ) 


In order to understand what is typical in the new way of putting 
the problem, we shall begin with the simplest case: the closed cylinder 

of infinite length which shows neither edges c = const., nor edges 
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s = const. Since for such a cylinder oblique waves can go around with- 
out interference, the appropriate expression for solution is easily 
guessed: With 


w = ^ a rL sin ^(ns + px) 


( 6 . 2 ) 


00 jr 

0 = XI An sin ij(ns + Px) 


where 2U is the circumference of the complete cylinder, the pair of 
equations (6.1) is transformed - with consideration of the expres- 
sion (3-1) and of the trigonometric identity 


2 cos sin ^(ns + px) = sin ^j((n+ 2r) s + Px ) + sin ^((n - 2r) s + Px) 


into 


+ n2 ) 2A n + ^ p2a n) 8ln ^ ns + ^ 


00 

r* 




n 

n 

^2r a n 

CQ 

H* 

P 

^(n + 2r)s + 

n=-°° 

r=l 


00 

r* 




r 

H 

^2r a n 

sin§ 

((n - 2r)s + 

n=-a> 

r=l 



y-L 

E 

_t 2 * 2 , 

(n2 _ L , 

„2\ 2 „ 1 , 

Y. L 

n=-oo\J- 

r 

- v c 

? U 2 

V* T “ ' °n ~ 1 
R o 


co r* 

H Y1 £ 2 r A n aln §(^ n + 2r ^ 3 + PX ) + 


n=-oo r=l 


00 r* 

H YY ^2r^n sin 5( (n - r)s + Px) = 0 
n=-co r=l 


( 6 . 3 ) 
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The equation system ( 6 . 3 ) is constructed exactly like the system (3 A) 

except for the fact that sin is now replaced by sin —(ns + Px), 

and P a by 2PnT. Nevertheless the further treatment is very 
different in the two cases . The summand px in the argument of the 
sine function causes the terms with positive and negative n no longer 
to be distinguished simply by the sign^ the expression (6.2) must there- 
fore contain all n values from -«> to 03 (whereas in the compression 
case the positive n values had been sufficient for obtaining a solu- 
tion that satisfied all conditions), and the equation system for deter- 
mining the a n and a_ n (/ -a n ) extends, therefore, to infinity "in 

both directions." That, moreover, the factor of the eigen-value 
is 2np instead of P 2 has the result that the minimum value for 1 
as a function of P lies in an entirely different p range than in the 
compression case. In order to perceive this last fact which is of 
utmost importance for the computational evaluation of the equation 
system ( 6 . 3 ) we shall here, exactly as in the compression case, briefly 
treat the cylinder of constant curvature at the outset. 

For £ 2r = 0 the equations (6.3) are each reduced to the first 

line. The requirement that the factor of sin £(ns + P x ) must disappear 

results in a system of equations which contain in every case only 
unknowns A n , a n with the same index. The buckling condition results 

precisely as it did in the compression case (compare equations (2 A), 
etc.) and reads with the abbreviations (2.5)* (2.5'), and 


= 8 — -A- = 4 1/3(1 - v2) 1 R 

CJ* /CD V ' * + 


E t 


according to (2.6) 


1 , n _ 

+ = ° 


(6.4) 


The condition = 0 yields 


0-J-+X , efl k 2) jnAn) g£ 

*n n n J d(3 2 « 


dP £ 




(6-5) 
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r - \/l2 fl - v2) u2 1/12(1 - v 2 ) 

Since for the complete cylinder /<d = — — — ^ = — — 

is a very great number, P must be either very large or very small in 
order to cause 


1 (n 2 + P 2 ) 2 + p3j/a^ 
n p y/oT (n 2 + p 2 ) 2 


(6.4') 


to assume reasonable amounts. With the assumptions 


P 2 « n 2 , k n 



there follow from (6. 5) for P 


and k n 


the simple formulas 



and therewith for M- 



*n */3 


(6.5' ) 


^ l) “ 






The critical shear load therefore attains its minimum for small P 
values. According to equation (6.2), p/n is the tangent of the angle 
formed by the node lines and the x-axis ; thus P « n signifies that 
long wave crests almost parallel to the axis develop. In case of longi- 
tudinal compression, on the other hand, wave fields with an aspect ratio 
of the order of magnitude unity are formed ) that is, the reciprocal 
action of stretching and bending, which is characteristic for the shell, 
takes full effect in case of compression, but not in shear; hence the 
critical shearing stress according to equation (6.6) is incomparably 
lower than the critical longitudinal-compressive stress according to 
equation (2. 8). 


"''■''For longitudinal compression there is ~ (t/R)^, for shear 

^ ( t/H) 3/2, for normal pressure (the wave crests run exactly axis- 
parallel) the critical hoop stress becomes ~ (t/R)^. 


c+ia 
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It is clear (and is confirmed by computational spot checks), that 
this behavior of the infinitely long shell under shear cannot change 
fundamentally, due to the variable curvature . In discussion of the 
system ( 6 . 3 ) we may, therefore, make use of the assumption 

(3 2 « n 2 

As a consequence, the quantities 



decrease very rapidly with increasing n, thus the secondary terms in 
the determinant of buckling for large n become negligibly small so 
that the discussion of the equation system resulting from equation ( 6 . 3 ) 
may be limited to a few terms around n = 0. Equation (6.6) is main- 
tained as an approximation result,* it permits the conclusion that the 
terms n = t2, as the terms with the smallest n, will play the decisive 
role; for, according to equation (6.4), n = 0 has no meaning, and 
n = 1 leads to an entirely different type of buckling (Greenhill 
buckling) with which we do not want to deal here and for the treatment 
of which the simplified equations (6.1) would not be sufficient. 12 

Thus we consider the system ( 6 . 3 ) for even n and limit ourselves, 
for reasons of clarity, to a three-term curvature expression (3*1) and 
to five equations. One then obtains by eliminating the A^, arranging 

and putting equal to zero the factors of 


sin ^y(ns + 0x) for n = -4, -2, 0, 2, 4 


a system of five 

equations. 

> the 1 

determinant of 

the coefficients of which 


(-4,-4) 

(■ 

- 4 ,- 2 ) 

( 

-4,o) 

( 

-4,2) 

( 

-4,4) 



(-2,-4) 

(■ 

- 2 ,- 2 ) 

( 

- 2 , 0 ) 

( 

-2,2) 

( 

-2,4) 


A = 

( 0,-4) 

( 

0 ,- 2 ) 

( 

0 , 0 ) 

( 

0,2) 

( 

0,4) 



( 2,-4) 

( 

2 ,- 2 ) 

( 

2 , 0 ) 

( 

2,2) 

( 

2,4) 



( 4,-4) 

( 

4 ,- 2 ) 

( 

4,0) 

( 

4,2) 

( 

4,4) 



-^Compare Fiugge: Statik und Dynamik der Schalen (Statics and 

dynamics of shells), p. 199 ff., Berlin 1935; furthermore, Kappus: 
Elastizitatstheorie grosser Verschiebungen, (Elasticity theory of large 
displacements), ZAMM'19, 1939* P- 351 ff. 
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must vanish. The coefficients (i,k) of the determinant are 
(°> 0) - t 0 + k o + *2% + *k\ 

= _l + k + 1 2 (* + k + 5 2 /k + k ' ) + a* 

(2,2) I kp 2 2 \ o 4 \ 2 6/ p 


(-4,. 

(4 


+ k 4 + ^( k 2 + k 6 


,-m] _ 1. 

,4) | k 4 

(-2,2) = (2,-2) = -(C 2 2 k 0 + 2kgC 4 ) 


( k 2 + H) + ^4 2 ( k o + k 8) 


+ V 

3 


(-4,4) = (4,-4) = -£ 4 2 k c 


-(-2,0) = -(0,-2) = (2,0) = (0,2) = C 2 (k 0 + k 2 ) + ^(kg + k^) 
-(-4,0) = -(0,-4) = (4,0) = (0,4) =^ 2 2 k 2 + 5 4 (k Q + k 4 ) 

(-4,2) = (2,-4) = (4,-2) = (-2,4) = -^ 2 C 4 (k o + k 2 ) 

(-4,-2) = (-2,-4) = (4,2) = (2,4) = ^ 2 (k 2 + k 4 ) + y 4 (k 0 + kg) 


Thus the determinant is built in a manner permitting an easy survey; 
by appropriate addition of the lines or columns it may readily be simpli- 
fied so that the very large term 

k = ^ 


appears only at the point (0,0). Since k , k,, k^ 

4 2 4 b 


are, compared to 


k o> small in the ratio E— , a four-series determinant remains which c^n 


n^ 


be further reduced if one also neglects everywhere k 4 , kg compared to 


k^ and makes sure that the large term 


_1_ = 2^6 
k 4 3 2 \/o d 


occurs only in the elements ( -4, -4) ( 4, 4) . The condition A = 0 is 
then reduced to the statement that the inner two-series determinant 
must vanish. In this manner a quadratic equation for p originates in 



28 


NACA TM 1302 


which the linear term also vanishes (the sign of n remains undeter- 
mined, as it must be), so that one has as buckling condition: 

^ + k 2 (l + S 4 - 2?/) 2 (6.7) 

This formula does not vary if a multiterm expression for l/R is 
used, because the £g, £g, in the decisive two-series determinant occur 

only in combination with k 4 , kg . . ., and thus must be eliminated 
corresponding to the other neglections leading to (6.7). 

The sign of £ 2 is insignificant (as it must be for the closed 

cylinder); the deviation from circular form causes a reduction in crit- 
ical load which can easily be calculated for the ovals of the type (4.8) 
or (4.9) that are of interest. First we must determine from 


( 6 . 7 ') 


the critical p -values by differentiation with respect to 3 or 1/3. 
The simple calculation yields for 1/3^ the quadratic equation: 

$ - [(* - r - 

3 ifs^ 1 ' £l *) 2 (l + Sit - 2?2 2 ) 2 ] " 0 (6.8) 

The two expressions 

(l - £ 4 ) 2 and (l + i k - 2£ 2 2 ) 2 

for the ovals of type (4.8), (4.9) in the £-range which is, according 
to figure 7> of interest, differ numerically only slightly so that it 
is permissible to replace, for the determination of the extreme, the 
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square of the geometric mean occurring in the last term of the equa- 
tion (6.8) by the square of the arithmetic mean; (6.8) then reads 

256 _ 2 (l - £h) 2 + (1 + Zlj. - 2C2 2 ) 2 _ 

o>P^ 2 


0#^ 

l 1 - ? 0 

2 + 

{ 1 + \- A i) 

f 

256 



2 



and has the positive root 


256 _ 1 „ t 1 - ^ 2 ) g + ( x + ^ - ge 2 g ) ; 


2\2 


o)P^ k 2 2 


= 3 


= 3^1 + (jj 2 - 2 ^ 2 2 ( 1 + + 2 ^ 2 ^) 

If we substitute this in equation (6.7)* the two radicands may be 
written in the form 

^(l + C^ 2 - 2^2 2 (l + + 2^) + 2(^4 - C 2 2 (l + £4) + $ 2 4 ) 


We have therefore 


H = 


+ ^ ■ 2y(1 + g + 2?2 T / t 1 ■ 

1/2 


k - s 2 2 ( 1 + sp + \ 2 

K 1 + i k 2 - 2^ 2 2 (l + Z k ) + 2^ 


(6.9) 


The last factor in (6.9) may be replaced by 1 for an oval of the type 
(4.9)- Then it can be written in the form 

<2VJ Z 


1 - ml !■ 


(l - hi 


1 + ^4 2 - 2£ 2 2 (l + £4) + 2^ 


which makes it clear that even in the extreme case 2^ = 0.6 X 0.9 
(compare fig. 7) the reduction caused by the second summand stays below 
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1 percent. Thus there results, with consideration of U = H Q n, for 
the critical shearing stress 


_ a*V a> 

T " ^ 


(l - v 2 ): 


^) 3/2 [ 1 + ^ + «!,) + V ]‘ 


(6.9') 


Aside from the factor & ♦ . • a i/k the formula ( 6 . 9 ') is built 
exactly like the formula for the circular cylinder derived by Flugge,- 
except that the numerical factors deviate slightly. According to 


Fliigge ' s theory the factor 


= O .236 replaces 


= 0 . 272 , thus 


j. u j. jr wnv ^ ^ xr~ . — — w • u. | 1 ^ 

3/2 3/3 

a numerical value smaller by 15 percent. This error in the 
formula ( 6 . 9 ') stems from the simplifications in the initial equa- 
tions ( 6 . 1 ) which were permitted only under the assumption that "many" 
waves developed in the hoop direction; however, the deviation is still 
within technically tolerable limits, in spite of the buckling form of 
extremely low wave number (n = 2 ). 

Thus the mean curvature l/R 0 ^ B ) according to (6.9' )> essentially 

decisive for the load capacity of the cylinder with oval directrix under 
shear load. One determines t/R 0 from the spar height 2b and the 

axis ratio e = — , with the aid of the e' -curve of figure 4 or 7 
a 

(according to the form of the directrix) on the basis of the 


relation — 
R~ 


2e' — J e' 

2 b 


is read off as a function of e by inter- 


polation of the variable £. Accordingly, (6.9') can be written in 
the form 

T / t\3/2,, , 


I = (±) 

E \2 b/ 


In figure 11 the factors ty(e), which - with consideration of Fliigge's 
correction - are given by 


*(0 = 


(1 - 0 . 09) 3/4 3/2 


— ( 2 e ') 3/2 1 + i k 2 - 2^ 2 2 (1 + + 2 ^^ 


^Schaleribuch (Book on shells), p. 20 6 . 
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are drawn (dot-dashed) for the ovals (4.5) and (4.9)* The two curves 
practically coincide. The expression ^(e) = 1.34c - 0.37^ - 0.25 
may serve as approximation formula for both ovals (and of course also 
for (4.8)); thus 



(6-9") 


7. THE BUCKLING CONDITIONS FOR THE CYLINDRICAL SEGMENT UNDER 
SHEAR LOAD; ITS SOLUTION IN FIRST APPROXIMATION 


More important for the application than the closed cylinder under 
shear load is the cylindrical segment of variable curvature supported 
along two edges s = constant under shear load (and shear-compression) . 
We think for instance of a half oval with the longitudinal edges 
supported according to the requirements of the conditions (2-9)- The 
expressions (3*2) and (6.2), which for the previous problems had led to 
finite-term equation systems, are not usable: One does not satisfy the 

differential equation, and the other cannot be adapted to the boundary 
conditions. In contrast, the problem leads to an infinite equation 
system in which every equation has an infinite number of terms, with the 
aid of the expressions^ 


w 



^ nx 

sin — + b^ cos 
i 11 



, ruts 

sin 

U 



, TtX 

sin — + B cos 
Z 



sin 


nrts 

U 


(7-1) 


which satisfy term by term the boundary conditions (2.9) • The cal- 
culation expenditure is nevertheless smaller than in the case of com- 
pression because the critical range is again the region of small 
(3 -values where the system converges extremely rapidly. For obtaining 


x Ifhese expressions are also used by Kromm In the report mentioned 
in footnote 5> Luftf . -Forschg. 1938; P- 517* Kromm gives a description 
of the physical importance of the expression (compare the formulation 
of (7-8) below in the summary) in the yearbook 1940, p. 834. 
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the equation system itself, one must employ considerations similar to 
those in section 3j it is particularly useful to introduce again the 
stipulation (3-6). If we limit ourselves, for reasons of clarity, 
preliminarily to the simplest curvature expression 


1 1 /, d.T\ s 

— = — 1 - 21 cos 


2ns' 


R R r 


U 


there result from the first equation (2.1) (compare (3*7)) 


TtX ICS 

by equating to zero the factor of sin P — sin n — 

U U 


An = ' (fJ 2 ! r?f + S ”-2)) 


(7-2) 


TCX 3TS 

and by equating to zero the factor of cos p — sin n — 


= - 


P^ 

EU^/ 


(f> 2 + 

.2 rn 

2 t . ' 

I n% 

H b n+2 + b n-2jj 


The second equation (2.1) where one can at first not factor sin 

U 

yields correspondingly 


1 


t2jt2 /o 2 


1 - V 2 U 2 


(p 2 + n 2 ) 2 - op 2 ) 


a n sm 


mis 

U 


R 


U 


— (l - 2£ cos 2it - ]P 2 5Z A n sin + 


U 


00 

y 

E 

z_ 

1 

1 - V 2 




(P 2 + n 2 ) : 


a P £ 


) b n sin - — (l - 2£ cos 2n — \p 2 ^ B sin - 
'J U E n \ U/ i n U 


2TpJI na n cos = 0 


(7.2') 


1 
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In order to make of these two equations - in which we visualize A n , 

B n inserted according to (7*2) - an equation system for determination 

of the a n , b n , we must replace the factors cos nTt i by sin nKS , by 

U U 

means of an additional Fourier expression. We first rename the summa- 
tion index in the last term 


2 T( 3 ^~ mb TO cos 

U 


The Fourier coefficients of the function 
the known manner from the relation 


cos 2 L_i are obtained in 
U 


■'U 


cosm^=V M 
U li- 


mits . nits , 
cos - 77 - sin - 7 - ds . 

U U sin £££ = it y . sin 321S. 

/ U sin 2 Eli ds U it V n 2 - m 2 U 

f o u 


Qm or n odd] 

The last term of the first equation (7*2') becomes therefore 


8 


o ? — c — , mn 

Jt n 2 . m 2 


m n 


sin 


nits 

U 


By exchanging the sequence of the two summations, introducing ( 7 . 2 ), 
putting the factors of sin £££ individually equal to zero and 

enlarging them by —— — , there results from the first equation ( 7 . 2 ') 

ybo 

' X ) &n + &n ( kn + ^( kn+2 + k n- 2 )) - ^ a n- 2 (k n + k n _ 2 ) - 

^ a n+2( k n + ^+2) + ^ 2k n-2 a n-k + ^K+2 a n+b + ~ “ ~~j=Yl ~2 ^2 = 

' It 3 a*V 01 m n c - ssf- 

(7.3) 

From the second equation (7*2') follows an identically constructed 
equation in which aj_ and b^ appear exchanged compared to ( 7 * 3 ), 
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and the sign of the last term is reversed. One now readily recognizes 
that one does not need both equation systems for determination of the 

critical value of T. 1 ^ For the second system is transformed into the 
first if one multiplies the first, third, fifth . . . equation by (-1) 
and puts 

^n = a n> ^m = ( - l) 


both systems have, therefore, the same eigen-value 


In order to clarify the further 
besides the abbreviations (2.5') and 

quantities s^ v by the equations 


considerations, we introduce 

M- = — additionally certain 
cr*/c I d 


■u - f ^ + d + » 2k i + - x ] 

S 22 * l + k 2 + - X 

s 33 = 5 + k 3 + ?2 ( k l + k ?) ‘ X 

s 13 - + + k 3 ] 

S = T 

15 ^ 3 

S 35 = P [ k 3 + ^ etc ’ 


(7- ft) 


15 The respective consideration also is that of A. Kromm, Luftf.- 
Forschg. 1938> p- 521. 
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Then the coefficient scheme of the equation system ( 7 . 3 ) is 


a l 

b 2 

a 3 

\ 

a 5 

b 6 

a 

7 

s 

11 

1.2 

s 

13 

1.4 
15 b 

S 15 

35 b 

0 . 

2.1 
T 1 1 

8 

22 

2.3 

s 

23 

21 

s 26 

2.7 

S 13 


s 33 

3^ 

7 

S 35 

3^ 

27 

S 37 ' 

iiin 

15 

3 2k 

7 

s 44 

9 b 

s 46 

4.7 
‘ 33 * 



(7.5) 


We investigate first the case of pure shear load (X. = 0). Pre- 
cisely as in the case of the complete cylinder, here also the lowest 
buckling load pertains to the buckling form with the "steep" waves, 
that is, to small (B -values^ as there, we need use, therefore, only a 
small initial portion (from the system (7>5)) and can make appreciable 
simplifications within the coefficients themselves. 

We obtain a first approximation for the critical values of ji 
from the determinant of buckling of the two first equations 



The pertaining wave length results from 


d (^I 2 ) 

d (3 2 ) 


-3( 8 ii 8 22 ) =0 


4 dP 


(7.6) 


(7.7) 


Because of 


(3 « 1 

the evaluation of equation (7*7) is very simple. From 

k i - ^ *2 - jg ^ \ = ^ e 2 /* 
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it follows first that the terms t^k^, £ 2 k^ cannot be of importance 

compared to their neighboring terms; if we, moreover, regard k, instead 

of p 2 , as the unknown to be determined, equation ( 7 . 7 ) is reduced to 
the simple problem of determining the minimum of the function 





fr + (i + ^St 1 + s 


as a function of kp. One finds that the second summand of the second 
parentheses can be omitted, since the differentiation of the first 

— (1 + £) 2 , which becomes zero for 


parentheses yields 


k, = P 2 y/m = 


1 + i 


(7-7') 


so that kp ^/256 in the neighborhood of the extreme is indeed vanish- 
ingly small. With (7«7')> and omission of the small terms, equa- 
tion ( 7 - 6 ) reads 


^ 2 = ;S e2 £ +(1 + 5)2k ^^-^ u + ?) 


'f3 2 \Jod 2 \f(E 


( 7 . 8 ) 


For the critical shearing stress T = ^1— p = . 

8 


tn |i 


1 - V^j U 2 


there results therefore 


T = 


3* 


-J-37l(?f 2 + ? - °' 55 E (ff /2 fi 7 * (7 - 9) 


2J2 (l - v 2 )' 

wherein, we can, of course, again write l/R 0 for jt/u. 
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8. CRITICAL SHEAR LOAD FOR THE HALF CYLINDER 
OF ARBITRARY CURVATURE 


The formula (7-9) has exactly the same construction as the 
formula (6.9') for the complete cylinder. However, it is not worth 
while to refine it or even to discuss it more closely because it 
represents, after all, among the formulas obtainable from the equation 
system (7-5) only the first approximation. Calculation of the second 
approximation does not offer any difficulties since one may now utilize 
the experiences made in determining the first. 

From the requirement that the determinant 


8 11 

3" 

S 13 

2 

S 22 

5 

3" 

S 13 


-b 

S 33 


has to vanish, the second approximation for the critic al value of 

p 

M- becomes 


( 


s ll s 33 " s 13 2 ) s 22 


’13 


II 4 36 8 

o S 33 + 25 S;L1 + s S !3 


- s i;L s 22 


8 h s 33 


We can put for the 


8l S - ] 1 l8 S-j 
1 + — + 3 




25b 33 5 S 


33 


( 8 . 1 ) 


3 ll “ + ( 1 + £) k i)> a i3 - - ^ P(l + 5)^, 


( 8 . 1 - ) 


3 22 


it J3_ _ 4jt|3 

4 1 ^ ’ 


3 = £ JL _ §ln 3_ 

33 4 k 3 " 4 k x 
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The correction terms in (8.1) thus become 


2 ¥(l + i)2 81 s n 


3 13 K 1 
81 


s u s 33 


1 + (l + ^) 2 k 1 2 ^ 25 


’33 


— [l + (1 + 0^ j. 


— — = - ^ (1 + « kl 


5 s 


33 


We can estimate their magnitude with equation ( T • T ' ) 


3 13 


1 


s ll s 33 


81 s n 18 s 13 2 


2t 


162(1 + £) 2 ' 25 s 33 5 S„ ~ 25 45(1 + 0 


33 


(8.2) 


One can see that for determining the locus |3 2 of the extreme value 

the numerator correction term may be eliminated; the denominator terms 

also may be cancelled, since the term 2/25 is a constant, and the 

remainder is so small that it is again of no importance. Thus — — = 0 

9 dP 

and for determining k-j_ and 0 , respectively, we retain the simple 

formula ( 7 • T ' ) : 


V - ± 

1 1 + g 


Thus equal signs appear in equation (8.2); by substituting equation (8.2) 
in equation (8.1), one obtains for the correction C to be made in the 
first approximation p-r 2 the expression 


1 - 


C = 


r 

162(1 + 


1 + TS? - 


2_ 

25 


2g 

45(1 + t) 


25 

27 


1 - 




162(1 + C) 2 _ 25 


l - 


10 

2^3 


1. 


1 + 1 


z* C- 
27 - 


(8.3) 


The quantities -12- 


5 


and 


ii 


affect Pjj ~ \fc in the 


243 1 + 5 162 (1 + o 2 

extreme case £ = 0.5 by less than 1 percent. The only essential 
difference between the two approximations p-j-j and Pj consists, 

therefore, in the numerical factor , /22. which causes reduction by about 

V 27 
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2 percent; thus a further approximation can not bring any further 
essential improvement.-*-^ If we now introduce t = p in 


II “ V 27 C 2* i I with 


C 2 - 


1 + 


£ 


i 2 — « 1 ♦ 1 « 


243 1 + i 324(1 + t ) 2 50 1 + i 


„ 3« . 

and h_ = -= -4 — — r-, we obtain 

1 /2 


T = 


for v = 0.3* thus 


W 6 (12(1 - E v 2 )) 3/ 2 ^”) V2 


? - °-533(f) 3/2 CjfTl 


(8.4) 


(8.4') 


For the special case £ = 0 this formula is transformed, as it must be, 
into Kromm's formula (valid for arbitrary u/R 0 ) : 


= 1,67E U|^ 


Ro 

If we introduce instead of U the spar height 2b, (8.4') reads 


(8.4") 


= 1.50(^(«) with \|r(e) = cp(t) = (e'U )) 3 / 2 C 2 (l + Z) 1 / 2 ( 8 . 5 ) 


wherein t is a function of the axis ratio which in the special case 
of the oval (4.5) is obtained from the curves of figure 4. It is drawn 
in a dashed line in figure 11. 

For an arbitrary oval represented by the general expression (4.1), 
the calculation is hardly different from that for the oval represented 

■^For the special case £ = 0 A. Kromm, Luftf . -Forschg. 1938, 
p. 525 > has proved by a check calculation that hjjj = 
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by (4.5)- The quantities in the more general case 5^, / 0 

are now, except for the factor ^ (3 , the quantities denoted by ( ik) 

(5.5) . The considerations made so far show that one may cancel a great 
deal; one obtains 



■13 - - ? 4 )(i - y*i 


( 8 . 6 ) 


and hence 



_w Sg - S 4 \ 

9 25 1 ' i6i\i - r 2 ) 

= 4 S 11 S 22 ^ l \ 10 ^2 - ^4 ' 

243 l - 5 2 


25 

% — S, , S 


12 


11 22 


1 k - k \ 

25 1 - 5 2 / 


(8.7) 


that is, 


T 

E 




1 £2 - C4 \ 

50 1 - 5 2 / 


(8.7') 


One can see that precisely as for the closed cylinder the sixth harmonic 
has no significance whatsoever, and that the difference compared to the 
special oval ^ ^ = 0 consists only in a modification of 

the unimportant last correction term. For the ellipse-like ovals we 
considered especially, the change means a very slight increase of the 
critical values, since according to equation (4.9) the factor of the 
correction term assumes in the limiting case 2^ = 1 the 
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value 


2- 1 = 1.2 


instead of 


= O .67 according to equation (4.5) 


1.75 1.5 

More important for the comparison between equations (8.7') and (8.4') is 
the fact that the factor 1 - for equal £-value is larger this time. 

For, according to equation (4.9), = -1.5C, and if one would simulate 


the oval still more to an ellipse (compare p. 16) , one would obtain 
coefficients still somewhat larger (1.6 for a curvature expression 

starting from sin® ££, etc.). 


The physical significance of this is that among the ovals of equal 
length U (and therewith of equal mean curvature 1/R q ) the buckling 

stiffness increases with the nonuniformity of the curvature variation. 

Of course this applies only if the spar supports the point of minimum 
curvature. If one places the support at the ends of the large axis 
(£2^0), it becomes more and more ineffective with increasing "ellip- 

ticity." We compare for instance the equations (6.9) and (8.7') with 
one another by substituting, according to equation (4.9), an< ^ 

in the extreme case 2£ = 1 the numerical factors 

0.291 X 0.72 = 0.21 for the unsupported cylinder 

and 0.533 \/o.25 X 0.99 = 0.26 for the supported cylinder 

result which now differ only slightly. If one intensifies the ellipticity 
by a curvature expression of still more terms (thus with still larger 
Cp values), the numerical values approach one another still more. 


If we again turn from the parameters appearing in equation (8.7') 
U and £ 


to — and 


e, we may write 


with 


T 



*U) = U'U)) 3/2 (i - y l/2 


1 _ ^2 - ^4 
50 1 - 


( 8 . 8 ) 

this time determined 


from figure 7- The result of the calculation is also plotted in 
figure 11 (dashed). Like _the curves pertaining to the closed 
cylinder, the two curves i|r(£) for the ovals (4.5) and (4.9) also lie 
very close together; that means that, as before in the compression case, 
the buckling load of the half oval supported at the "long sides" is 
almost not at all dependent on the exact form of the oval, but only on 

the two parameter values — and — = e . Since the relations (8.4), 

t a 
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k2 

(8.7'), etc., as the comparison of the exact Flugge formula with 
formula ( 6 . 9 '), contain (due to the simplifications made in the initial 
equations) possibilities of error of an order of magnitude of 15 percent, 
it is sufficient to replace the two T curves found here (as well as 
all others possible) by the curve drawn solidly in figure 11. This 
approximation curve has the simple equation 

♦ Q (e) = 3-le - 1.3e 2 - 0.3 (8-9) 

To repeat this formula, as the result of the calculations of the 
two last sections, once more: the simple relation 

^ 3//2 (3.ie - l-3e 2 - 0.3) (8.9') 

is valid for the critical shear load of the infinitely long statically 
supported half oval. It is so simple that it is not worth while to 
represent it by a diagram. According to its derivation its validity 
range is limited to the region 1^6^. l/3; however, most probably not 

too large errors will arise if it is applied for the region e ^ 

to e * i. 

5 


9. THE WING NOSE UNDER COMBINED SHEAR 
AND COMPRESSIVE LOAD 


The extraordinary diversity of the buckling form makes it impossible 
for compression and shear buckling to "mix"; in case of an appropriate 
variation in load the one buckling form changes suddenly into the other, 
that is, the curves for the critical values in a a-, T_coordinate 
system show a break. (Compare fig. 12a.) The calculation which must 
be performed according to the methods indicated in section 5 shows that 
the critical compressive load within very wide limits is changed not at 
all by Bhear, whereas, as we shall now demonstrate, presence of com- 
pression or tension influences the critical shear load. 

We start from the equation system (7*5)> thus consider first again 
the special oval £4 = = . . . = 0. The first approximation has 

again the form (7*6); however, this time the quantities s must be 

IK 

understood to signify the complete expressions (7-^) (that is, with 
X / 0) . Since the presence of a small compressive force is not able to 
basically change the critical behavior of the cylinder, we may in the 
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Vcjd region which is of interest for the wing nose again substitute 

r 2 

k., = — j- y/ 7 B and obtain (compare also ( 7 . 8 )) 
n 4 




16 


k 16 k 


= 9« 2 /l 


+ (1 + 0 2 ki - X 1 - 


Xkl' 

IT; 


(9.1) 


For determining the critical 3 values from 

d (^I 2 ) _ .nr. d (^I 2 ) = 0 


d3‘ 


= 


dk„ 


Xk-, 

the small term may be omitted. Therewith there results for the 


critical value of k^ the value (7-7 1 ) independent of X: 

k i = p 2 ^ = rfr ; 


( 9 . 2 ) 


\ij thus becomes: 


p 2 = 2d[i + jl x+ _£_ 

1 2 >/S 7 \ 16 32(1 + £), 


(9-3) 


For the plotting of p against X the tangential direction at the 
point X = 0 is of foremost importance; there results from the simplifi- 
cation of (9*3) which is valid for small X values: 



9X \ 
32(1 +5)/ 


(9-3') 


Exactly as in the case of pure shear load a corresponding formula 
results, without considerably larger calculation expenditure, also for 
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the second approximation pertaining to the system (7.5) 
again from the formula (8.1): 


One starts 


II 


9 

— s s 
4 11 22 


1 - 


13 


B ll s 33 


fill . fil 

25 833 5 833 


( 9 . 4 ) 


where one has to substitute for s 


ik 


1 


11 = ^ F" + (1 + ^ k i “ x > s l 3 = - % p(l + 


22 




\ R 1 


16 


33 4 81 


( 9 - 5 ) 


The term B l^ f B ±± 9 ^ ^' tie numera tor may be omitted, as we have 

seen above (p. 38). Then one obtains equation ( 9 . 4 ) in a readily 
surveyable form if one multiplies throughout in the numerator and the 

Xk, 

denominator by 1 - — = 

81 


“11 2 - 


fe * (1 * Q\ - x)(^- - jg)(i - 
(' ■ 17) + 5 ( 1 + f 1 + »V - ^)-|u * «h 2 


( 9 . 4 ') 


The position of the minimum of this expression lies again close to 

. _ 1 

k i _ r+i 


If one substitutes this value, one obtains 


2 9 it 2 25,- 1, v 

^TT = ~ -7= — (1 + t) 

11 2 27 


l 1 ~ 2ii | nX 1 ' iSii IttK 1 ' 81(1 + l)) 


1 - 


106 X 


10 5 


21B7 1 + £ 2F3 1 + 5 


(9.6) 
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For the tangents at the point X = 0 one obtains by series development 
(compare also equation (8.4)) 




50 1 + s; 


1 - 0.263 


1 + 1 


( 9 . 6 -) 


or written in t and o, with t(X = 0) denoted as t 


t(\) = T 


“ — 

1 n . A 


~ cr2b 0.263\/3(l - V 2 ) 

1 + 5 

" T o 

Et € '(i + i) 


(9.6”) 


The formulas (9*6) to (9.6 ") are valid for the special oval (4-5). The 
corresponding formulas for the general oval (4.1) differ from them by 
the fact that instead of 1 + £ one has to write 1 - £ 2 and that the 

very small correction term ^ — or — £ — , respectively, 

243 1+? 50 1 + i r _ £ 

contains, instead of the second fraction, the quantity — _. In 

1 - k 

determining it had been found that these formal differences are 

almost completely compensated at transition to the two parameters 2b/t 
and b/a, so that for t q the expression (8.9') which is independent 

of the oval form resulted. With still better approximation one may 
replace the factor representing the influence of the longitudinal force 


0.263\/3(l - v 2 ) = 0.433 = 0.433 

e'(l - £ 2 ) e'(l - y t*( 6 ) 


"by a "mean" function independent of the oval form. In figure 11 the 
function \|/*( e) is drawn for instance for the two ovals (4.5) and (4.9). 
The deviations are indeed so slight that - in view of the simplified 
initial equations - it would he perfectly nonsensical to take them into 
consideration. 


In place of (9*6") we obtain therefore the formula, valid for 
arbitrary ovals, for the critical shear in the presence of small 
compressive or tensile forces (a is as compression counted positive): 


T 


0.433/ g gb\ 
t*U)\E t / 


(9-8) 
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On the basis of the formulas (9-8) (for small X.) or, respectively, 
(9-6) (for larger X values) one may sketch a diagram from which result 
the admissible critical pairs of values a, T. In figure 12 a repre- 
sentation was chosen which saves use of separate diagrams for deter- 
mination of the critical loads t q , ct q (everytime in the absence of 

the other) . This is possible because of two peculiarities of the curves . 
The dependence of the curve branches starting from the r-axis ( shear 
with some compression) on the parameter 2b/t is such that one may 
factor this parameter, as equation (8-9') shows, by selection of the 
easily calculated abscissa 



The remaining one-parameter curve family is then numbered according to 
the parameter e = iL. The curve branches starting from the a-axis 

(compression with some shear) cannot be transformed, in a similarly 
simple manner, into a one-parameter family dependent only on e, by 
selection of a suitable ordinate. Instead of this, however, these 
curves are independent of T , that is, they are horizontal straight 
lines which need not be explicitly drawn into the diagram; both param- 
eters on which - if the ordinate is selected - these "curve" 

E 2b 

branches depend can now be indicated by a simple method: one need only 

mark the pertaining 2b/t value at the intersection of the horizontal 
straight line and the other branch of the a-, T-curves numbered 
according to e. Connecting, in addition, the points of equal 2b/t 
values, one has on the whole a diagram with two one-parameter curve 
families from which one may read, without additional auxiliary diagrams, 
the relations between the four quantities a , t, 2b/t, b/a which 
characterize the critical state. 

For use of the diagram (which of course is drawn only for the 
region 50-1000 of the parameter 2b/t that is of practical interest) 
the following method results: b/a is known, t is selected, 2b/t 

determined. The point of intersection of the curve pertaining to the 
respective b/a with the abscissa axis represents the pertaining 
critical t q value in the absence of compression (formula (8.9')). 

For smaller t compression is permissible; it is read off from the 
curves indexed b/a vertically above T. If the point comes to lie 
in the region covered by the dot -dashed curves, one must consider 
whether or not compression alone produces instability. This is the 

case when the intersection of the curve indexed b/a with the normal 

2b 

comes to lie above the dot-dashed --curve, since the piercing 
points of the dot-dashed with the b/a curves indicate the critical 
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compression values in the absence of shear. Only if these "critical 
curves" are not reached, shear is permissible. The method is best made 
clear by an example 


- = 0.60 




1 . — (— = 1 . 0 : a compression — = 0 .l 6 — is permissible 


E 


I fa\3/2 , k: 

E \t / 


2 b 


a compression — = 0.29 — is permissible 
E 2b 


since the dot-dashed curve —■ = 400 lies higher up; 


3 -i T 


2b \ 3/2 


0.8o. One would read off — = 0.57 tut this 

E 2b 


value is not admissible, since at £ = 0.34 75 -- lies the critical value 

jj 2b 

for compression alone. 


If the permissible compressive and shear loads prescribed by the 
problem lie higher (lower) than the critical ones found from the diagram, 
the reading must be repeated with a larger (smaller) t value. 


10. VALIDITY LIMITS - UNSYMMETRICAL OVALS 
UNDER BENDING - SUMMARY 


The report investigates the stability of the cylindrical shell of 
infinite length and variable curvature simply supported at the longi- 
tudinal edges in case of loading by longitudinal compression and shear. 
Starting from the cylinder equations ( 2 . 1 ) which are reduced to the 
"essential" terms, the formulas are developed first in such a general 
form that the calculation may be performed for a shell of arbitrary 
curvature (representable by a Fourier polynomial) and of arbitrary opening 
angle and then be numerically evaluated for the special ovals (4.5)* 

(4.8), and (4.9) with the opening angle n (wing nose). It is found 
that the shell under compression shows an entirely different behavior 
from that of the circular-cylindrical shell (the oval shell of equal 
mean curvature buckles very much earlier), that, however, in case of 
shear loading - where almost axis-parallel wave crests originate - the 
variability of the curvature has only little effect. Performance of 
the calculation is merely more troublesome in case of compression than 
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in case of shear, because the equation system to be solved in the param- 
eter region of importance for the wing nose converges there very much 
more slowly. On the other hand, the result in case of compression 
(fig. 10) has a very much wider range of application - since the waves 
in longitudinal direction are short, the results are, with good approxi- 
mation, valid also for the cylinder of finite length. For shear, in 
contrast, presupposition of great cylinder length is essential - if, 
due to the transverse reinforcements, the oblique waves can no longer 
go around without interference, a constraint results which sometimes 
increases the critical load to a multiple. It is difficult to tell 
by how much one may, however, assume that it amounts to less than for 
the circular* cylinder because for shortened longitudinal -wave length 
the variability of curvature (insignificant for the cylinder of infinite 
length) must have a reducing effect, as in case of compression. Besides, 
for the breaking load one may count on a considerable reserve also in 
case of compression; noticeable buckling originates at first only in 
the slightly curved part, without much affecting the strongly curved 
part, so that a considerable amount of load may be shifted to the latter 
before the entire shell collapses. (The shell of strongly variable 
curvature behaves similarly to a plate stiffened by longitudinal supports 
above the buckling limit.) 

The diagrams and formulas have been calculated under the assump- 
tion that the ovals are symmetrical to the large axis, that they show 
the opening angle cp Q = it, and that they undergo compression or tor- 
sional load. However, they may be used also when the ovals are unsym- 
metrical, have an opening angle different from it, and undergo bending 
with or without transverse force. The investigations of section 5 
(compare fig. 9) have demonstrated that the two halves of the oval show 
almost no mutual influence across the strongly curved nose. One cal- 
culates, therefore, an unsymmetrical oval under compression by mirroring 
the flatter one of the two halves (as the more endangered one) at the 
large axis and taking the dimensions of this substitute oval as a basis 
for the calculation. If the opening angle is cp Q / it, one provides an 

equivalent half oval by slightly modifying the more endangered half of 
the "apexes" so that there the tangent is vertical to the end tangent 
at the flattest place and mirrors the thus originating quarter oval; 
this is, as an approximation, permissible, because the region of smallest 
curvature is of foremost importance for the critical behavior so that 
a slight change in the region of greatest curvature (which practically 
remains unbuckled) cannot essentially affect that behavior. If, 
finally, the wing nose does undergo, not pure compression (bending of 
the wing about the vertical axis), but bending about the large axis. 


"^Compare the circular-cylinder results of Kromm, Jahrbuch 1940 
der deutschen Luftfahrtforschung, p. 1832. 
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one mirrors the compressed half and calculates the thus originating oval 
as if the maxi mum bending-compressive stress were uniformly distributed. 
With this method, one is on the safe side without calculating very 
unfavorably: for even for the circular cylinder the critical bending 

stress is not very much higher than the critical compressive stress,!® 
and for the oval the difference is a great deal slighter still, due to 
the width of the pressure zone (which is moreover the least curved) . 

Transverse forces will, in general, not affect the stability of 
the wing nose. Vertical forces will not do so because the part trans- 
ferred by the wing nose is smaller and, moreover, reaches its maximum 
value in the strongly curved portion, and horizontal ones will not, 
because the stresses in the nose must change the sign and thus can 
generally not be large. Therefore only the torsion remains as essential 
shearing load. In case of combination of torsion and bending about the 
horizontal axis, the diagram 12 serves for the half under compression. 
For the half under tension (which sometimes, nevertheless, can be the 
more endangered one because it usually is the flatter half) we obtain 
from (8.9) and (9*5) the approximation formula 


T = 


'1 + 


0.^33 fa 2b\\ 

**(«) U t )) 


in which cr, as tension, has to be inserted as positive. The quan- 
tity f*(e) is taken from figure 11 - for rough calculations one may 
put f* = 1. 

We add a compilation of the abbreviations and of the end formulas 
important for practical use. (Compare also fig. 2.) 

Abbreviations: TT 

3=y (2-3) 


\f03 = 


llhl - V 2 ) 

* 2 v' 


e t£ 
1 - v 2 3 u 2 


( 2 . 5 ) 


_1_ = 1_ (p 2 + n 2 ) 2 

k n p2 


x = 2 / 3(1 - v 2 ) £ 5a, 


, o ^ : v2) 1 ^ 


(2.5') 




l8 Flugge, Ing . -Arch. 3, 1932, p. 501 ff. 
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For the half cylinder U = R Q n and thus 


^ = Vi g ( 1 - v2 I 2 = 2/3(1 - v2 ) ^ = ^3(1 - v 2 )^— 

n t t t 


Formulations: 

for the curvature: 


— = — (l + 2^2' cos — - + 2^ cos + . . . + 2^ r * cos 


R R, 


2ns 


4ns 


2r*ns 


U 


U 


U 

(4.1) 


for the unknowns w and 0: 
(a) for compression 


00 

v . nx , 
w = / sm — sin 


u . nns 
U 


j . nx . nns 

0 = 2_ An srn t sin -fr 

A L U , 

(h) for shear (closed shell) 

00 " 

v a n sin ^ |3x + 118 ) 

-00 

0 =X A n sin ^-(px + ns) 
“00 

(c) for shear (shell segment) 


( 2 . 2 ) 


( 6 . 2 ) 


w = 


n=l,3,5 


o nx _. nns , 

a n sin — sm + 


b„ cos — - sin 


nx nns 


— — -r ✓ ouo — — 

u n .2A6 n 1 


u 


a X a nx , nns , X ^ ^ roc . nits 

0 = Z_ An 8in T sin ~lT + ^ , B n cos T Bln IT 

n=l, 3>5 n=2, 4, 6 


> (T.l) 


The critical loads are represented as functions of the two param- 
eters 2b/t and b/a (compare fig. 1) in figures 10 (pure compression) 
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and 12 (shear and compression). For pure shear approximation formulas 
can be given: 

(a) complete cylinder 


i = (^) 3/2 ( 1 - 34e ■ 0<3Te2 ■ °- 25 ) (6 * 9 " j 

(b) half cylinder 



Translation by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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(c) Ellipse. 




Figure 7.- Expression (4.9). For conversion of the parameter £ 

1 ^ . ^ 
to € - — (axis ratio) and e' = = ~ (ratio between small 

axis and mean curvature radius), R^, R 2 are the apex curvature 
radii. 
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Figure 9.- Compression buckling profile plotted over half the developed 
width U/2 of a half shell; shown in the secondary figure are the 
amplitudes of the seven harmonics of which the buckling pattern is 
composed (curvature expression (4.5)). 
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Figure 10.- Critical compressive stress a as a function of the 
parameters 2b /t and b/a. 
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